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Abstract 

The interaction of a quantum deformed oscillator with the environment is 
studied deriving a master equation whose form strongly depends on the type of 
\1J \ deformation. 

O : PACS number(s): 03.70. +k, 02.50.-r 
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QlhI Quantum groups [|^ introduced as a mathematical description of deformed Lie 

algebras have given the possibility to generalize the notion of creation and annihiliation 
C \ operators of the usual oscillator and to introduce g-oscillators H, H . Soon after, several 
attempts have been made to give a physical interpretation of the g-oscillators. In Ref. 
qh| Q they appear in connection with the relativistic oscillator model, while in Ref. 
g-oscillators have been used in the generalized Jaynes-Cummings model. 

Instead, the nature of g-oscillators of the electromagnetic field is clarified by the 
^ ■ nonlinearity of the field vibrations 0. This g-nonlinearity reflects the exponential 
H i growth of the frequency of vibrations with the amplitude. 

This observation suggests that there might exist other types of nonlinearities for 
which the frequency of oscillation varies with the amplitude via a generic function 
/; this leads to the concept of /-deformed oscillators 0. Then, the generation of 
coherent states for a class of /-deformed oscillators enter in the real possibilities of 
trapped systems 0. 

By virtue of possible further physical realizations of /-oscillators, it could be inter- 
esting to study the environmental effects on such oscillators. 

The aim of this paper is to describe the behavior of a nonlinear oscillator, like the 
/-oscillator, plunged in a bath modeled by an assembly of harmonic oscillators. A 
master equation approach is developed to consider several types of reservoir, restricted 
however to the case of small damping. 

The free Hamiltonian of an /-deformed oscillator may be taken, in the case of 
unitary frequency and natural units, as 

h, = ]^{a^a + aa^) , (1) 



where the operators A and A^ result as a distortion of the usual annihiliation and 
creation operators a and a^, 

A = af{n) = f(n + l)a, n = a)a^ (2) 

A^ = f{n)a^ = a^f{n + 1) . (3) 

Here, / is intended as an operator- valued function of the number operator (moreover it 
is assumed Hermitian). In general, it can be made dependent on continous parameters, 
in such a way that, for given particular values, the usual algebra is reconstructed. This 
is the case of g- deformations |0 . As a consequence of the noncanonical transformations 
(|^) and (^), the commutation relation is not preserved, 

A,A^]={n + l)f{n + l)-nf{n). (4) 

Of course the usual algebra is restored whenever / — ;> 1. Furthermore, the Eq. (0) can 
be rewritten as 

no = l[nfin) + in + l)fin + l)] , (5) 

from which results more clear its nonlinear character. 

Now, the usual model of reservoir consisting in a large number of harmonic os- 
cillators is employed by means of the Hamiltonian Hb = Y^j^jbjbj and the linear 



interaction between system and bath |jTO| 



n^nt = art + ^tr . (6) 

The bath operator F is given by F = J2j i^jbj, with Kj the frequency dependent coupling 
constants. In write down Eq. @, we have assumed that the frequency of bath modes 
most strongly interacting with the system is very large compared to the strength of the 
interaction (measured in frequency units). This assumption is also known as Rotating 
Wave Approximation [ p!0|] . 

We imagine that the interaction (^) is switched on at the time t = and the bath 
at this time have no correlations with the oscillator. The initial joint density operator 
u'(O) can thus be taken as the product of some unspecified operator p(0) and a given 
operator p_b(0) for the bath. It obeys the Liouville equation 

W{t) = -% [Hq + nB+ Hint , w{t)] = Lw{t) , (7) 

where the superoperator L should be intended as the sum Lq + Lb + L^nt- The reduced 
density operator of the oscillator will be 

p(t) = TrB{e^*u;(0)} = f/(t)p(0). (8) 

Assuming the inverse of the operator U{t) to exist, we may write p(0) = U{t)^^p(t), 
and thus the formal equation of motion 

p(t) = f/(t)f/(t)-VW^£(t)p(t). (9) 

Restricting ourselves to the weak damping case, C{t) can be constructed perturbatively 
by an expansion in powers of the interaction part Lint of the Liouvillian superoperator 



L. Hence, neglecting the third and higher order terms, the reduced equation of motion 
is obtained 

p{t) = Lop{t)+J^ dt' I (r(t')rt)[a(-t')p(t),at] + (^t^(t'))[a^ p(t)a(-t')] 

+ (r(t')r)[atp(t),at(-t')] + (rtrt(t')) [ap(t),a(-t')] 

+ H.c.j, (10) 



where the expectation values are taken on the bath and represent its correlation func- 
tions. 

To go on, it is also necessary the explicit time dependence of the operators a and 
a^; from Eq. (^), it is 

a{t) = exp [—iQ{n)t] a (11) 

with 



n{n) = -\{n + 2)f{n + 2) - nf{n) 



(12) 



an operator- valued frequency constant in time. 

In Eq. (10), by assuming a short correlation time of the bath operators (shorter 
than any typical system periods), the extreme of integration t can be replaced by oo. 
Furthermore, it is possible to invoke the fluctuation-dissipation theorem in the following 
form 



oo 



dt{V{t)T^)e 

rft(rtr(t)) 

dt{VV{t))e 



Ant 



int 



x{n)N{n) 

x(fi)M(fi) 



where 



xm = g{n)K\n) , 



(13) 
(14) 
(15) 

(16) 



is the bath response function determined by the density function g of the bath modes 
in the continuum limit [|l^]. Hence, the master equation (p!0D becomes 



p{t) = -t[no,pit)] 

+ ^p.v. f dn' x{^')l {i + N{n')) [{n-n')-^ap{t),a^ 



+ N{Q') a\p{t){Q-Q'y^a + M{Q') a^ p{t),a\Q - Q') 



v\-i 



+ M^{Q')\ap{t),{Q-n')-^a 
1 ' 



H.c. 



+ 



1 + 



iV(fi))x(fi)ap(t),at] + [at,p(t)iV(fi)x(fi)a 



+ 



aV(t),a^M(n)x(fi)l + \ap{t),M\^)x{^)(^ + H.c. 



(17) 



where P.V. denotes the Cauchy principal value and the following relation has been used 



/ 



2 27r 



dn'{n~n')-^j^{n') 



with J-' a generic bath correlation function and JF its Fourier transform. 

Now, the Markov approximation x(^) = 27 |l^, can be introduced in the second 
term with curly brackets in Eq. (0), while it is inadmissible in the terms involving 
the principal value frequency integrals since they become divergent. The latter are the 
analogous of the Stark and Lamb shifts in the common theory of damping [|l^]. To 
simplify the treatment, however, these terms are disregarded from now on, due to their 
small contribution. 

A suitable basis to represent the master equation (p!7[) in c-number form is obviously 
the number states basis, then one gets 



Pm,n\J') 



mf{m) + (m + l)f{m + 1) - nf{n) - (n + l)f{n + 1) 



7< m 



l + N{Vt{m-l)) +nl + N{Vt{n-l)) 



Pm,n\J') 



+ (m + l)iV(fi(m)) + (n + l)N{n{n)) p^,„(t) 



+ 7V(m + l)(n + 1) 2 + N{n{m)) + N{Sl{n)) p„+i,„+i(t) 



+ '-)\/nin 



N{n{m - 1)) + N{n{n - 1)) p„_i,„_i(t) 



+ 7 Jm(n + 1) M{n{n)) + M{Q{m - 1)) p^_i,„+i(t) 



+ 7^(m + l)n[M*(n(m)) + M*(fi(n - l))Jp„+i,n-i(t) 

- 'y^Jm{m - l)M{Q{m - l))p„_2,n(^) 

- -f^n{n-l)M*{n{n - l))p™,„_2(t) 

- 7y'(m + l)(m + 2)M*(fi(m))p„+2,„(t) 

- 7^^(72 + l){n + 2)M{n{n))p,^^n+2{t) , 



(19) 



with pm,n{t) = {^\p{t)\n). 

It is worth noting how the system deformation strongly affects the structure of the 
master equation ( [T9| ) also in its irreversible terms through Q. Instead, usually the 
nonlinearities are included in the master equation only through the term Lq, so that 
the damping terms remain unaffcted by the nonlinearities. This seems incorrect, unless 
M{n) and iV(fi) go to zero. 

A further simplification is made by assuming to have only a thermal bath, so that 



M(fi)=0; N{n) 



1 



expipn) 



(20) 



where f3 = l/ksT; ks is the Boltzmann constant and T the temperature. In such a 
case the stationary solution of the Eq. ([T9|) can be easily found 



Pm,n (00) = 5m,nZ CXp 



P 



-^{nf\n) + {n+l)f\n + l)) 



(21) 



where Z is a normalization constant. The solution ( pTD may result far from a trivial 
Boltzmann distribution when a deformation is present; for example, in the case of 
a (small) g-type nonlinearity it immediately gives the deformed Planck distribution 
derived in Ref. [|ll|]. 

Finally, it is to remark that the derived master equation is not of the Lindblad form 
||12||, and this inevitably leads to situations in which unphysical density operator will 



arise (one illuminating example could be given by the harmonious deformation [|T3 
for which f{n) = \j ^Jn). Nevertheless, this is not a reason to abandon this equation, 
it only means that one has to be careful in applying it, evaluating the range of validity 
for the kind of considered deformation. 

Summarizing, a quite general master equation describing the damped dynamics 
of a deformed oscillator has been proposed, showing how its structure reflects the 
nonlinear character of the oscillator. The specific solution strongly depends on the 
type of introduced deformation. On the other hand, it could be also searched by 
transforming the master equation in a partial differential equation for a quasiprobability 
distribution [Q, exploring possible connection among deformation parameter (s) and 



ordering parameter(s) [|14| 

Acknowledgments 

I express my gratitude to V. I. Man'ko for introducing me to the subject of deformed 
algebras. 



References 

[1] V. G. Drinfeld, Proc. Int. Congress of Mathematicians, 1, 798 (1986); M. Jimbo, 
Lett. Math. Phys. 11, 247 (1986). 

[2] L. C. Biedenharn, J. Phys. A 22, L873 (1989). 

[3] A. J. Macfarlane, J. Phys. A 22, 4581 (1989). 

[4] R. M. Mir-Kasimov, in Lecture Notes in Physics, Vol. 382. Proc. 18th Colloquium 
on Group Theoretical Methods in Physics, (Springer, Berlin, 1991), p. 215. 

[5] M. Chaichian, D. Elhnas and P. Kuhsh, Phys. Rev. Lett. 65, 980 (1990). 

[6] V. I. Man'ko, G. Marmo, S. Solimeno and F. Zaccaria, Int. J. Mod. Phys. A 8, 
3577 (1993). 

[7] V. I. Man'ko, G. Marmo, F. Zaccaria and E. C. G. Sudarshan, Proc. 4th Wigner 
Symposium, N. Atakishiyev, T. Seligman and K. B. Wolf Eds. (World Scientific, 
Singapore, 1996), p.421; Physica Scripta, 55, 528 (1997). 

[8] R. L. de Matos Filho and W. Vogel, Phys. Rev. A, 54, 4560 (1996). 

[9] A. P. Ploychronakos, Mod. Phys. Lett. A 5, 2335 (1990). 

[10] see e.g., W. Louisell, Quantum Statistical Properties of Radiation, (Wiley, New 
York, 1989). 



[11] V. I. Man'ko, G. Marmo, S. Solimeno and F. Zaccaria, Phys. Lett. A 176, 173 
(1993). 

[12] G. Lindblad, Commun. Math. Phys. 48, 119 (1976). 

[13] E. C. G. Sudarshan, Int. J. Theor. Phys. 32, 1069 (1993). 

[14] K. E. Cahill and R. J. Glauber, Phys. Rev. 177, 1882 (1969). 



